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We investigated the optical properties of epitaxial n- and p-type ZnO films grown on lattice-matched 
ScAlMgCU substrates. Chemical doping yielded a severe inhomogeneity in a statistical distribution of involved 
charged impurities. Two approaches are adopted to treat the inhomogeneity effects; Monte Carlo simulation 
technique for the n-doped films, and the fluctuation theory for the p-ZnO. The broadening of PL band of n-ZnO 
was significantly larger than predicted by theoretical results in which the linewidths of each individual emis- 
sions have been determined mainly from the concentration fluctuation of donor-type dopants by the simulation. 
Moreover, the rather asymmetrical line shape was observed. To explain these features, a vibronic model was 
developed accounting for contributions from a series of phonon replicas. In case of p-type ZnO:N, analysis 
of excitation-intensity dependence of the peak shift of donor-acceptor luminescence with a fluctuation model 
has also proven the importance of the inhomogeneity effect of charged impurity distribution, as in the case of 
ZnO:Ga. We extracted the inhomogeneity in the sample and acceptor activation energy prepared under the var- 
ious growth conditions. It is shown that the theoretical results are in good agreement with the experimental 5-K 
time-resolved luminescence for the systems in a fluctuation field. Finally, localized-state distributions have been 
studied in N-doped ZnO thin films by means of transient photocurrent measurement. 



I. INTRODUCTION 

Among widegap II- VI semiconductors, ZnO seems to be one of the most promising materials for optoelectronic applications. 
This is due to its stable excitons, having a large binding energy of 59 meV HHH]. It is also desired to more deeply understand 
the underlying doping physics yl |4|] for the achievement of a compound-semiconductor-quality p-n junction [ 5 ] . In our previous 
work, we reported the observation of intense near-band-edge photoluminescence (PL) for ZnO:Ga at room temperature (RT) [ 6] . 
It can be easily imagined that a high doping level tends to yield a severe inhomogeneity in a statistical distribution of involved 
charged impurities. In this work, such an effect is characterized with the Monte Carlo simulation technique |7|]. The phonon 
interactions are also discussed with a 'vibronic' model [8|]. 

Donor-acceptor pair (DAP) photoluminescence (PL) bands have been widely used to characterize acceptor-doped semiconduc- 
tors. Especially, in wide-gap semiconductors, it has been difficult to obtain /?-type doped materials. However, the time-resolved 
PL has not hardly been studied M lldl for ZnO:N. For the case of ZnO, a large amount of nitrogen well in excess of 10 19 cm -3 
has to be doped to achieve /?-type conductivity even using one of the best growth techniques that have ever reported [a repeated- 
temperature modulation (RTM) technique] [5]. It has become apparent that the "standard" theory formulated by Thomas, 
Hopfield, and Augustyniak ifTTll (hereafter, it is called a "standard" THA model) is totally incapable of explaining results of the 
DAP PL for heavily doped materials probably because this does not take the effects of inhomogeneity into account. In this work, 
to implement the effect of the fluctuation field inherent to the doped samples, we successfully grew /?-type-dopants-incorporated 
films, which are suitable for the experimental verification of the relevant "fluctuation" theory. 

II. EXPERIMENTAL PROCEDURES 

We have grown gallium- and nitrogen-doped ZnO films grown by laser-molecular-beam epitaxy on lattice-matched 
ScAlMg04 substrates. The detailed methods of fabrication have been given elsewhere [5]. A very complicated growth se- 
quence has been adopted for the nitrogen doping. We measured the photoluminescence (PL) of the doped ZnO layers in the 
band-edge region. The experimental procedures [6] were identical to those adopted in our previous study and their details were 
reported elsewhere 11211 . Pulsed excitation was provided for the time-resolved experiments with the frequency doubled beam of 
a picosecond mode-locked Ti: sapphire laser. 
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III. RESULTS AND DISCUSSION 

A. Optical properties of n-type doped ZnO 

FigureQJa) shows PL spectra of ZnO:Ga at four different doping levels (dashed lines). The shape of the PL band is asymmet- 
rical. The linewidths of the PL increase from 154 to 195 meV with an increase in Ga concentration (n&j). The "theoretical" 
broadening for the PL was evaluated in terms of potential fluctuations caused by the random distribution of donor impurities 
(cf. Fig. 3 of our previous paper J6[]). These values are listed in Table I. However, the experimentally observed broadening are 
greater than the calculated widths (<7i). 




PHOTON ENERGY (eV) 

FIG. 1: (a) Room-temperature photoluminescence spectra (dashed curves) of n-type ZnO doped with different Ga concentrations. Also shown 
are the results of the fit to the data using a 'vibronic' model (solid lines), (b) The lowermost PL curve of frame (a) with individual contribution 
of the emission lines (dash-dotted lines). 



TABLE I: Sample specification and main characteristics for the four «-ZnO samples. Ga concentration (nc a ), theoretical broadening ct\, 
broadening evaluated from the Monte Carlo simulation 02, zero-phonon peak energy E, and the deduced Huang-Rhys factor S. 
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It is thought that the theoretical widths are exceedingly underestimated because the model might be too much simplified. We 
here adopt another approach properly taking the microscopic fluctuation of donor concentration into account; the thermalization 
redistribution model developed by Zimmermann et al 11311 . According to this model, there is a relationship between the Stokes 
shift and the full-width at half maximum (FWHM) of PL. One can estimate the energy scale of the band potential profile 
fluctuation from the Stokes shift and predict consistently the FWHM. Figure|2]shows a Stokes shift of the luminescence plotted 
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FIG. 2: The PL Stokes shifts (closed squares) plotted against the Ga concentration. The measurements were performed at RT. 



against the gallium concentration. Obviously, the shift energy increases with incorporation of the donor dopants. The energy 
scale of the band potential profile fluctuation can be evaluated from the Stokes shift and predict consistently the FWHM. A Stokes 
shift of the luminescence is plotted against the gallium concentration as shown in Fig. [2] At sufficiently high temperatures, the 
above-mentioned relationship is indicated as Estokes = —<J 2 /kT. The nomenclatures (fee and T) take their conventional 
meanings. For example, the Stokes shift for the sample with the nc a of 8.0 x 10 18 cm -3 is 22 meV, leading to a of 24 meV. 
It also yields the temperature independent value of 2ay / (ln4) ~ 46 meV in the FWHM, which can be also supported by the 
Monte Carlo simulation. 

Miller- Abraham's rate for phonon-assisted carrier tunneling between the initial and final states i and j with the energies of -E 1 , 
and Ej was adopted to simulate the hopping events of photo-created carriers: 

2r i3 {Ej- E 3 + \Ei- Ej\) \ 

a 2kT J ' ' 

Here is the distance between the localized states, a is the decay length of the wave function, and Vq is the attempt-escape 
frequency. Hopping was simulated over a randomly generated set of localized states with the sheet density of N. Dispersion of 
the localization energies was in accordance with a Gaussian distribution, 

g(e) aexp(-(e-£ ) 2 /2c7 2 ). (2) 

With the peak positioned at the mean resonant energy Eq and the dispersion parameter (the energy scale of the band potential 
profile fluctuation) a. The radiative recombination event with the probability t^ 1 terminates the hopping process. The emission 
energy at the time of recombination is scored to the luminescence spectrum. By implementing the thermal broadening effect, 
the overall FWHM of the PL band is given by: 

FWHM (ct 2 ) w [(2,r v /(ln4)) 2 + {1.8k B T) 2 } 1/2 . (3) 

Figure [3] represents, with a solid line, the simulated result obtained for the following values of parameters: Na 3 = 0.1, 
toVq = 10 4 , and a ~ 24 meV. The lineshape of resulting spectrum was independent of the choices of the Na 3 and tqUq 
parameters. The simulated curve could be well approximated by the Gaussian function, which is not surprising if considering 
the measurement temperature (300 K). The FWHMs (02) of the simulated PL spectra for the samples at four different doping 
levels are complied in Table I. 

The simulated broadening (02) is smaller than that of experiments as in the case of the simpler theoretical consideration. This 
difference is probably due to the contribution of the phonon replicas. In addition, the simulated line-shape is rather symmetrical, 
which is in great contrast to the asymmetrical shape observed in the experimental spectra. We think that the zero-phonon peak 
and its phonon replicas are not spectrally resolved well, leading to the larger experimental broadening. It is inferred that the 
coupling with longitudinal-optical (LO) phonon is relatively strong. ZnO may have coupling constant significantly stronger than 
that of GaN or of ZnSe Q. 
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FIG. 3: Simulated PL spectrum under stationary excitation. The Ga concentration is shown in the figure. 



Only the terms of longitudinal optical (LO) phonons are taken into account [ 15], whose interaction is known to be the strongest. 
In such a case, the transition energy of the replicas will be hu> = E'zpl — ^lo^lo where £'lo=72 meV in ZnO, -Ezpl is the 
energy of the zero-phonon peak and t^lo is an integer. A Gaussian function with width o\ (cf. Table I) modeled each individual 
emission line in the spectrum. The probability of a given phonon emission is proportional to (5^0° /f?LoO> where Slo is the 
Huang-Rhys factor for the LO phonons and ?7lo the number of phonons emitted in a transition. The solid lines in Figs.QIa) and 
(b) were obtained by summing all of the lines using an appropriate set of parameters. Table I summarizes the donor-impurity 
concentrations (riGa), zero-phonon peak energies (E), and S-factors. The respective contributions whose peak positions are 
shown by five dash-dotted curves in Fig. [Jib), with equidistant peaks by the E^o of 72 meV. We can say from this fit that the 
asymmetrical and broad PL band can be explained by the contribution of LO phonon replicas. It should be noted that, even in the 
undoped ZnO, the zero-phonon emission becomes weaker than LO-phonon-assisted annihilation processes at elevated sample 
temperatures. 



B. Optical properties ofp-type doped ZnO 

Thin ZnO films are known to be grown more or less under residual «-type. The trial of the p-type doping often results in an 
observation of the radiative recombination from donor-acceptor pairs (DAPs). Monte-Carlo simulation [7] cannot be used for 
the analysis of such a DAP system because the photo-carriers in these structures cannot be regarded as a single particle. The 
experimental PL features cannot be explained with the simple THA model in the case of our N-doped ZnO samples. Symbols in 
Fig- IS a) show a set of 5-K PL decay curves taken at three different emission energies, while the corresponding time-integrated 
PL is given in Fig. |4jb). We first remark that the decay time constants for all the energies are similar with respect to each other 
and the decay curves are non-exponential. This similarity is difficult to be predicted in the framework of THA theory. In this 
model, using the pair distance (R) between donors and acceptors, the emission energy (hv) is given byd: 

hv = hv x + e 2 /eR = {E G -E A - E D ) + e 2 /eR, (4) 

where the other symbols take their conventional meanings. Equation gives a unique relation between hv and R. The radiative 
recombination kinetics of DAPs is predominantly determined by the pair distance, R. It is well-known that the decay time is 
strongly dependent on R, and hence on hv. This strong "dispersion" effect is in apparent contradiction with our observed data. 

According to the treatment of Kuskovsky and his coworkers [ 16] in which the fluctuation effects are accounted for, an equation 
describing the time-dependent PL intensity is now quite different from that in its counterpart [11]. By averaging over all possible 
realization of the fluctuation potential, transient intensity /^(t) at an energy E is: 
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FIG. 4: (a): Comparison of the experimental time-resolved PL on ZnO:N at three different emission energies with a fluctuation theory. Open 
circles are the experimental data, while the solid curves are the theoretical fits with parameters of £ = 13, rj = 0.036, n = 0.011, and 
W^max = 3 x 10 10 s -1 . The PL intensities were normalized and each curve was vertically shifted for clarity. The measurement temperature 
is 5 K. Secondary-ion mass spectroscopy (SIMS) was used for the determinations of the total N concentrations ([N] ~ 6 x 10 18 cm -3 ). Its 
conductivity is not yet p-type. Also shown in (b) is the corresponding time-integrated PL. 
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The zero energy as hvoo, defining an effective emission energy corresponds to E = hv — hfoo- We have introduced a di- 
mensionless time t — W max t and energy E = eRsE/e 2 , as well as parameters £ = 2R s /Rb, T) = e 2 / (AkT g eR s ), and 
n = tt(Na — Nr>)R%, which have allowed us to have all integrals as dimensionless factor. Here, Na an d No are the acceptor 
and donor concentrations, R s is the screen radius, and T g is the frozen temperature. With the use of Eq. it is possible to show 
that the decay is non-exponential, relatively close to stretched-exponential, or even slower, particularly at very low temperatures. 

The solid curves of Fig. Ufa) refer to the results of theoretical fit at different hv according to Eq. (0. For the calculation, all 
the parameters except E were kept fixed (£ = 13, r\ = 0.036, n = 0.011, and W max = 3 x 10 10 s _1 ). These values are quite 
reasonable, assuming Rb — 23 A (which corresponds to the donor depth of 30 meV) and R s = 150 A. As seen from the figures, 
the theory results in reasonably satisfactory fits for all hv's. The fitting procedure is sensitive to the values of the net acceptor 
concentration (Na — No) and of W max . 

Although the experimental verification for this theory has been performed in N-doped ZnSe, it should be emphasized that 
only samples influenced by the very strong fluctuation fields have been analyzed [17, 18]. In regard with this theory, wide 
applicability on the magnitude of fluctuations has not been yet verified, which is the main focus of the present work. This theory 
introduces a concept of screening radius (R s ) which is intimately related to the degree of fluctuations (smaller R s corresponds 
to the severer fluctuations). 

We discuss the time-integrated PL [cf. Fig. HJb)] from the viewpoint of excitation-dependent shift to quantify R s of our 
sample. The PL is dominated by the D S AP recombination (D s means a shallower donor) at 3.235 eV. The D S AP band exhibits 
no detectable shift with a change in excitation, suggesting very large radius: R s > 100 A. We confirmed that the use of 
R s = 150 A gives the well-reproduced excitation-intensity dependence of hv. Here, we used a value of 3.228 eV 11911 as 
hVoo = Eq — Ea — E D . E G = 3.428 eV, E A = 170 meV and E D = 30 meV JH. On the other hand, the similar analysis 
for ZnSe:N yielded in: R s ~ 20 A, indicating that our sample is appropriate to verify whether this theory can also be used to 
analyze the system with significantly larger R s . 

One of the reasons for the unexpectedly short decay compared to those in other semiconductors JH] is the difference in R s , 
because the decay becomes slower with increasing fluctuation. The other one is presumably related to the W max , which is one 
order of magnitude larger than that for ZnSe, and may be inherent to ZnO Xiong and coworkers 111 Oil have also reported the 
DAP PL lifetime of a ZnO:N bulk crystal shorter than many other semiconductors. The probability W max is represented by the 
following equation: 



128 e 2 iihv 00 E p 



3 moh 2 c 3 



oa_ 



(6) 



where /i is the index of refraction, mo is the electron mass, c is the velocity of light, E p is a band parameter that has been calcu- 
lated by Lawaetz using a k p analysis, and the result for ZnO is 28 eV, and the other symbols take their conventional meanings. 
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The calculation using the effective mass theory for effective Bohr radii of donors and acceptors evaluated the theoretical W max 
to be: W^ ax = 4 x 10 9 s" 1 , which is one order of magnitude smaller than the experimental value. It should be, however, noted 
that this theoretical evaluation sensitively depends on the choice of ZnO material parameters adopted in the calculation of the 
Bohr radii of donors and acceptors. We would like to say that the experimental value of E p is likely to be greater than 28 eV, 
the contention of which is based on the comparison of the values with other semiconductors and the strong excitonic effects of 
ZnO. Note, E p is related to the matrix element between the conduction and the valence bands. 



C. Photoconductivity spectroscopy of N-doped ZnO 



The ZnO:N films show a complicated transient behavior of the photoconductance in reaching their respective dark or illumi- 
nated equilibrium conduction state. Figure [2 a) shows a typical current buildup and decay illuminated with sub-band-gap light 
of 384 nm for ^20 minutes. The photoexcitation as well as the decay of photoexcited carriers are very slow. 

After turning off the light, the photo-induced conductivity is observed to prevail for very long times. It is evident that 
the decay is not essentially exponential for ZnO:N; a rather common feature for persistent photo-current (PPC) experiments, 
implying interaction between the photo-carriers and localized states during the decay processes. If the PPC is coming from the 
AX-center-like defects, the rising behavior at the edge in PPC excitation spectrum must take the form of (hut — Tiui^) 1 - 5 where 
hujth denotes the threshold energy. 




Time (min) Energy E ' (eV) 



FIG. 5: (a) Rise and decay transients of the current in ZnO:N at 293 K. The dashed line stands for a saturated current, (b) Density of states 
obtained from the PPC transient using a high-resolution LTM (open circles) in which the dark current was subtracted. The continuous line 
corresponds to the result of fit to these circles. 



To support this, we calculate the spectral distribution of the trap density-of-states using the multiple trapping model. Resulting 
rate-equations are given as follows; 

dn (t) dn,j(t) n(t) 

— = -l^— t r +n ° m > (7) 



= ujin{t) -jini{t), (8) 

where n is free carrier density, rij is carrier density captured at the i-th localized state, r is radiative recombination time for free 
carriers, no is photoexcited carrier density, uii(oc g(Ei) x AE) is a capture rate of carrier, ji is a releasing rate, Ei is an energy 
of i-th localized state, and g is density of state. 

By solving the Laplace-transformed rate equations analytically, one obtains relationship of g{E' ) with the Laplace transform 
of the PPC transient I{t, s): 

9(E'o)= mWU-^Wf, (9) 
yV 0) <jvk B T [7( s )]3 



E' Q = k B T ln( — ), 



(10) 
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where it is a capturing cross section, v is thermal velocity, s is a Laplace variable, v is a detuning frequency, Jit) = tl(t), 
K(t) = t 2 I(t), and E' Q is the trap energy below the mobility edge. The other symbols used in the equations have the usual 
meaning. 

The LTM thus gives a view of the density-of-state distribution g, as was shown in Fig. 12b). The spectral lineshape follows 
an exponentially decaying function. The former component was fitted using an exponential function of the exp(— ftw/fceTo) 
type, where To means the characteristic temperature. The other symbols have the usual meanings. The best fit to our data is 
obtained for To = 387 K. Comparing with the results on a polycrystalline ZnO films B22I1 . deduced spectral distribution of g in 
our case is completely different from that of the polycrystal. In the latter case, g had a peak-like distribution with a maximum 
near the one-third of the energy gap. This can be attributed to the large difference in the density of grain boundaries [23]. The 
characteristic temperature T deduced from Fig.Qb) is 387 K. 



IV. SUMMARY 

Gallium and nitrogen were successfully applied as sources of donor and acceptor impurities in the laser molecular-beam 
epitaxy growth of n-type ZnO, respectively. Their optical properties have been investigated in terms of their effect of inhomoge- 
neous dopant distribution and their phonon interaction. We also discussed the PL shifting behavior in ZnO:N and compared them 
with the controls grown by the conventional method. We have explained how the degree of dopant distribution inhomogeneity is 
extracted from the results of optical characterization. The activation energy of acceptors has been evaluated by proper modeling 
to be ~ 170 meV, which is independent of nitrogen concentration. We have observed transient behavior of photoexcited carriers 
in semi-insulating N-doped ZnO using photocurrent spectroscopy. The obtained persistent photocurrent shows decay times pre- 
vailing over hours. We obtain the spectral distribution of density-of-state distribution of the localized states. The lineshape of 
the spectrum suggests an existence of localized exciton states. The exponential fit [exp(— ^cj/fceT))] deduced a characteristic 
temperature T of 387 K. 
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